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Divergent  cusped-field  thrusters  are  electric  propulsion  devices  characterized  by  a  di¬ 
vergent,  dielectric  wall  channel  and  a  magnetic  topology  dominated  by  two  magnetic  cusps 
and  a  magnetic  bottle  near  the  anode.  Initial  testing  results  suggest  that  this  technol¬ 
ogy  may  lend  itself  to  miniaturization  while  providing  improved  lifetime  performance  over 
miniature  Hall  thruster  designs.  This  paper  presents  results  of  detailed  probing  of  the  di¬ 
vergent  cusped-field  thruster  plume  using  a  faraday  cup  and  a  retarding  potential  analyzer 
to  measure  the  ion  current  density  and  energy  distributions  at  varied  operating  conditions. 
These  results  provide  a  more  complete  experimental  characterization  of  the  plume  than 
was  previously  available  by  measuring  plume  attributes  at  wider  angles  with  increased 
spatial  resolution.  Plume  data  analyses  show  regions  with  the  greatest  current  density 
contain  the  most  energetically  uniform  ion  populations.  However,  the  retarding  potential 
analyzer  results  indicate  that  highly  energetic  ions  are  expelled  at  large  angles  relative 
to  the  thruster  axis  of  symmetry,  though  the  current  densities  at  these  angles  are  much 
lower.  Plume  measurements  made  while  varying  the  cathode  operating  condition  are  also 
presented.  The  operating  conditions  of  the  divergent  cusped-field  thruster  and  cathode  are 
found  to  be  strongly  coupled,  motivating  a  future  study  that  systematically  examines  the 
effect  of  cathode  location  and  operating  condition  on  the  performance  of  the  thruster.  The 
mode  transitions  that  occur  while  operating  the  divergent  cusped-field  thruster  are  also  ex¬ 
amined,  providing  a  quantitative  comparison  between  the  two.  Additionally,  time  resolved 
anode  current  measurements  are  presented  for  several  operating  conditions.  These  data 
reveal  a  distinct  difference  between  the  low  and  high  current  modes,  demonstrating  that 
the  anode  current  is  much  less  oscillatory  in  the  low  current  mode  than  the  high  current 
mode. 


Nomenclature 

Va  =  Anode  potential  with  respect  to  ground 
ma  -  Anode  mass  flow  rate 
Ik  =  Keeper  current 
mc  =  Cathode  mass  flow  rate 
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Figure  1.  DCFT  visualizations. 


I.  Introduction 

Stationary  plasma  type  (SPT)  Hall  thrusters  have  been  operated  in  earth  orbit  since  19721  with  typical 
efficiencies  near  50%  and  a  specific  impulse  of  1600  s.2  However,  the  operational  lifetime  of  SPT  Hall  thrusters 
is  often  limited  by  erosion  of  dielectric  material  used  to  provide  insulation  from  the  plasma.3  Researchers  at 
the  Massachusetts  Institute  of  Technology  (MIT)  Space  Propulsion  Laboratory  (SPL)  began  investigating 
divergent  cusped- field  (DCF)  thrusters4  6  to  determine  their  promise  as  a  longer  lifetime  alternative  to 
existing  low  power  Hall  thrusters.  Other  cusped- field  thrusters  include  the  High  Efficiency  Multi-stage 
Plasma  (HEMP)  thrusters  patented  by  THALES  Electron  Devices7  9  and  the  Cylindrical  Hall  Thruster 
(CHT)  developed  at  the  Princeton  University  Electric  Propulsion  &  Plasma  Dynamics  Laboratory.10  A 
detailed  comparison  of  the  DCFT  design  to  the  HEMP,  CHT,  and  other  cusped-field  devices  is  provided  by 
Courtney  et  al.6 

A.  Divergent  Cusped-Field  Thruster  (DCFT) 

DCF  thrusters  are  electric  propulsion  devices  characterized  by  a  divergent,  dielectric  wall  channel  and  a 
magnetic  topology  dominated  by  two  magnetic  cusps  and  a  magnetic  bottle  near  the  anode.  The  magnetic 
circuit  consists  of  three  annular  samarium  cobalt  (SmCo)  permanent  magnets  arranged  such  that  like  poles 
face  one  another,  and  a  core  at  the  base  of  the  thruster  composed  of  1018  carbon  steel  (see  Figures  1(a) 
and  1(b)).  Xenon  gas  is  injected  5  mm  upstream  from  the  downstream  face  of  the  anode  through  a  porous, 
stainless  steel  disc  embedded  in  the  magnetic  core.  In  this  way,  the  neutral  injection  and  anode  locations 
are  decoupled  from  one  another  in  the  DCFT  design.  For  complete  details  of  the  DCFT  and  existing  data, 
please  consult  previous  studies. 4-6,11-13 

Existing  experimental  data4-6  catalogues  DCFT  performance  (anode  specific  impulse,  thrust,  anode 
thrust  efficiency)  at  operating  conditions  ranging  between  100  W  to  500  W,  with  a  nominal  operating  condi¬ 
tion  corresponding  to  a  xenon  flow  rate  of  8.5  seem  and  an  applied  anode  voltage  of  550  V.  At  this  operating 
point,  the  DCFT  consumes  242  W  of  power,  and  produces  13.4  mN  of  thrust  while  operating  with  an  an¬ 
ode  thrust  efficiency  of  44.5%.  Although  basic  DCFT  operating  principles  are  thought  to  include  azimuthal 
electron  currents  within  regions  characterized  by  radial  magnetic  fields,  the  competing  roles  that  strong  mag¬ 
netic  bottles  and  internal  electric  fields  play  in  determining  the  location  and  efficiency  of  ionization  (whether 
in  the  upstream  cusp,  downstream  cusp,  or  both)  have  yet  to  be  explicitly  measured  in  experiments  or 
rigourously  explained  based  on  theoretical  arguments  or  simulation  results.  Demonstrating  the  effectiveness 
of  the  magnetic  topology  and  divergent  chamber  geometry  as  erosion  barriers  remains  unachieved  as  well. 
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The  present  investigation  provides  detailed  probing  results  of  the  DCFT  plume  using  a  Faraday  Cup  (FC) 
and  a  Retarding  Potential  Analyzer  (RPA)  to  measure  the  ion  current  density  and  energy  distributions  at 
different  flow  rates  and  applied  voltages.  FC  and  RPA  results  are  presented  for  several  operating  conditions, 
including  cases  where  the  cathode  flow  and  power  were  altered  to  study  the  affect  of  these  changes  on  visible 
and  measured  plume  structure.  Anode  current  oscillation  measurements  are  also  given. 

II.  Plume  Measurements 

The  following  experiments  were  performed  in  Chamber  1  at  the  Edwards  Air  Force  Base  (AFB)  Air  Force 
Research  Laboratory  (AFRL)  Space  Propulsion  Branch.  Chamber  1  is  a  cylindrical  non-magnetic  stainless 
steel  vacuum  chamber  with  a  2.4  m  diameter  and  a  4.1  m  length.  Pumping  is  provided  by  two  liquid 
nitrogen  baffled,  1.2  m  flanged  gaseous  helium  two  stage  cryogenic  (15  K)  vacuum  pumps  with  a  measured 
pumping  speed  on  both  xenon  and  krypton  of  48,500  L/s.  An  automated  2D  (R-0)  stage  was  used  to  move 
Faraday  and  RPA  probes  within  the  plume  of  an  operating  DCF  thruster  to  measure  the  characteristics  of 
the  plasma  being  ejected.  The  arm  mounted  on  the  rotary  stage  is  a  stepper  motor  drive  linear  translational 
stage  capable  of  traversing  up  to  80  cm  from  the  center  of  rotation  of  the  rotary  stage.  This  combination  of 
rotary  and  linear  motion  allows  for  the  continuous  sweeping  of  various  plasma  interrogation  probes  in  the 
R-0  plane.  Further  details  of  the  probes  and  experimental  facilities  are  provided  by  Nakles  et  al.1A  Figures 
3(a)  and  3(b)  show  the  thruster  firing  at  the  Edwards  AFB  AFRL,  and  Figure  2  illustrates  the  coordinate 
system  used  to  present  the  probe  dat.  Thermocouples  located  on  the  thruster  exterior  near  the  anode 
and  outer  case  enclosing  the  permanent  magnets  indicated  that  the  temperature  did  not  exceed  200°  C  at 
those  locations  during  all  tests  performed  for  this  investigation.  Measured  chamber  pressure  did  not  exceed 
3.5  /iTorr. 


R=  0  cm  I 


-0 


0  =  0 


Figure  2.  Top  view  of  the  2D  (R-0)  coordinate  system.  The  0  =  0  line  is  aligned  with  the  DCFT  axis  of 
symmetry.  The  cathode  is  above  the  probing  plane,  as  seen  is  Figures  3(a)  and  3(b). 


(a)  DCFT  operating  in  the  high  current,  or  diffuse, 
mode. 


(b)  DCFT  operating  in  the  low  current,  or  colum- 
nated,  mode. 


Figure  3.  A  visual  comparison  of  the  DCFT  plume  in  two  distinct  modes,  operating  at  the  Edwards  AFB 
AFRL. 
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Figure  4.  RPA  mounted  and  aligned  to  measure  ion  energy  distributions  within  the  DCFT  plume. 


A.  RPA  Probing 

This  section  presents  RPA  results  for  the  different  operating  conditions,  presented  in  Figures  5(a)  through 
9(f).  The  RPA  data  are  angularly  resolved  from  ±90°  or  ±120°  off  axis  in  5°  intervals  in  both  directions, 
taken  at  two  different  radii  from  the  thruster  exit  plane  at  each  condition. 


4  of  15 


The  31th  International  Electric  Propulsion  Conference,  Ann  Arbor,  Michigan,  USA,  September  20-24,  2009 
Distribution  A:  Approved  for  public  release;  distribution  unlimited. 


Retarding  Potential  (V) 


Retarding  Potential  (V) 


(d) 


(e) 


(f) 


Retarding  Potential  (V) 
(§) 


Retarding  Potential  (V) 

(j) 


Retarding  Potential  (V) 

(h) 


Retarding  Potential  (V) 

(k) 


Retarding  Potential  (V) 

(m) 


Retarding  Potential  (V) 

(n) 


Retarding  Potential  (V) 

(i) 


Retarding  Potential  (V) 
(1) 


Retarding  Potential  (V) 

(o) 


(8) 


Figure  5.  Ion  energy  distributions  plotted  at  varying  positions  with  respect  to  the  DCFT  exit  plane.  Operating 
condition:  Va  =  300  V,  ma  =  8.5  seem,  I*.  =  0.25  A,  mc  =  1.0  seem.  The  units  of  all  radial  distances  are  given 
in  centimeters.  r  f  ir 
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Figure  6.  Ion  energy  distributions  plotted  at  varying  positions  with  respect  to  the  DCFT  exit  plane.  Operating 
condition:  Va  =  450  V,  rha  =  8.5  seem,  I*.  =  0.25  A,  rhc  =  1.0  seem.  For  angles  outside  of  db  90°  at  this  operating 
condition,  see  Figures  8(a)-8(f).  The  units  of  all  rac|j^jf  cjiptances  are  given  in  centimeters. 
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Figure  7.  Ion  energy  distributions  plotted  at  varying  positions  with  respect  to  the  DCFT  exit  plane.  Operating 
condition:  Va  =  450  V,  rha  =  8.5  seem,  I&  =  0.25  A,  mc  =  1.0  seem.  For  angles  of  magnitude  greater  than  90° 
at  this  operating  condition,  see  Figures  9(a)-9(f).  Tl^fUji^its  of  all  radial  distances  are  given  in  centimeters. 
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Figure  8.  Ion  energy  distributions  plotted  at  varying  positions  with  respect  to  the  DCFT  exit  plane.  Operating 
condition:  Va  =  450  V,  rna  =  8.5  seem,  I&  =  0.25  A,  mc  =  1.0  seem.  For  angles  outside  of  d=  90°  at  this  operating 
condition,  see  Figures  6(a)-6(s).  The  units  of  all  radial  distances  are  given  in  centimeters. 


Figure  9.  Ion  energy  distributions  plotted  at  varying  positions  with  respect  to  the  DCFT  exit  plane.  Operating 
condition:  Va  =  450  V,  rha  =  8.5  seem,  I&  =  0.5  A,  rac  =  1.0  seem.  For  angles  of  magnitude  less  than  90°  at 
this  operating  condition,  see  Figures  7(a)-7(s).  The  units  of  all  radial  distances  are  given  in  centimeters. 
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Figure  10.  Ion  energy  distributions  plotted  at  varying  positions  with  respect  to  the  DCFT  exit  plane.  Operating 
at  Va  =  451  V,  iri a  =  4.3  seem,  I&  =  0.8  A,  rhc  =  2.0  seem. 
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B. 


Current  Density 


Figure  11.  Faraday  Cup  mounted  and  aligned  to  measure  current  density  within  the  DCFT  plume. 

This  section  contains  the  faraday  cup  results  for  the  different  operating  conditions,  plotted  as  the  current 
per  steradian.  The  data  are  given  in  Figures  12  through  14.  The  FC  data  are  angularly  resolved  from  ±120° 
off  axis  in  1°  intervals  in  both  directions,  taken  at  two  different  radii  from  the  thruster  exit  plane  at  each 
condition. 


-100  -50  0  50  100 

Angle  (Degrees) 


Figure  12.  Ion  current  density  distributions  plotted  at  varying  positions  with  respect  to  the  DCFT  exit  plane. 
Operating  condition:  Va  =  300  V,  rha  =  8.5  seem,  I&  =  0.25  A,  iinc  =  1.0  seem. 


10  of  15 


The  31th  International  Electric  Propulsion  Conference,  Ann  Arbor,  Michigan,  USA,  September  20-24,  2009 
Distribution  A:  Approved  for  public  release;  distribution  unlimited. 


Figure  13.  Ion  current  density  distributions  plotted  at  varying  positions  with  respect  to  the  DCFT  exit  plane. 
Operating  condition:  Va  =  450  V,  rha  =  8.5  seem,  I&  =  0.25  A,  rnc  —  1.0  seem. 


Figure  14.  Ion  current  density  distributions  plotted  at  varying  positions  with  respect  to  the  DCFT  exit  plane. 
Operating  condition:  Va  =  450  V,  rha  —  8.5  seem,  I&  =  0.5  A,  mc  =  1.0  seem. 


Figure  15.  Ion  current  density  distributions  plotted  at  varying  positions  with  respect  to  the  DCFT  exit  plane. 
Operating  condition:  Va  =  475.7  V,  rha  =  5.8  seem,  1^  =  1.48  A,  mc  =  1.0  seem. 
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III.  Current  Oscillation  Measurements 


This  section  presents  the  time  resolved  anode  current  as  recorded  with  an  oscilloscope  and  processed 
with  a  signal  analyzer.  They  demonstrate  a  conspicuous  difference  between  the  anode  current  behavior  in 
the  high  and  low  current  modes.  The  anode  currents  shown  in  Figures  16(a)  through  16(f),  representing  the 
DCFT  operating  in  the  high  current  mode,  are  characterized  by  high  magnitude,  low  frequency  oscillations. 
For  all  high  current  operating  conditions,  the  lowest  harmonic  of  the  frequency  spectrum  is  between  3.5  kHz 
and  4.0  kHz.  The  flight  time  of  a  neutral  Xenon  particle  with  a  temperature  of  1000  K  is  approximately 
0.2  ms.  The  frequency  corresponding  to  the  inverse  of  this  time  is  5  kHz.  The  anode  currents  shown  in 
Figures  17(a)  through  17(f),  representing  the  DCFT  operating  in  the  low  current  mode,  do  not  exhibit  the 
same  oscillatory  behavior  as  seen  in  the  high  current  mode.  The  large  oscillations  present  in  the  high  current 
mode  are  reduced  by  several  orders  of  magnitude.  Previously,4, 5  the  low  current  mode  was  found  to  be 
the  more  efficient  of  the  two  modes.  The  lack  of  strong  plasma  oscillations  in  the  low  current  mode  further 
suggest  that  it  is  the  preferred  mode  of  operation. 
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(a)  Time  resolved  anode  current.  Operating  condi¬ 
tion:  Va  =  300  V,  rha  =  8.5  seem,  I&  =  0.25  A, 
rhc  =  1.0  seem. 


Time  (sec) 

(c)  Time  resolved  anode  current.  Operating  condi¬ 
tion:  Va  =  450  V,  rha  =  8.5  seem,  I&  =  0.25  A, 
rhc  =  1.0  seem. 


Frequency  (Hz) 


(b)  Anode  current  spectrum.  Operating  condi¬ 
tion:  Va  =  300  V,  rha  =  8.5  seem,  1^  =  0.25  A, 
rhc  =  1.0  seem. 
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(d)  Anode  current  spectrum.  Operating  condi¬ 
tion:  Va  =  450  V,  rha  =  8.5  seem,  1^  =  0.25  A, 
rhc  =  1.0  seem. 


Time  (sec) 

(e)  Time  resolved  anode  current.  Operating  con¬ 
dition:  Va  =  450  V,  rha  =  8.5  seem,  I =  0.5  A, 
rhc  =  1.0  seem. 


O.OE+O  2.0E+4  4.0E+4  6.0E+4  8.0E+4  1.0E+5 

Frequency  (Hz) 


(f)  Anode  current  spectrum.  Operating  condi¬ 
tion:  Va  =  450  V,  rha  =  8.5  seem,  I&  =  0.5  A, 
rhc  =  1.0  seem. 


Figure  16.  Time  resolved  anode  current  and  spectrum  for  the  high  current  mode. 
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(a)  Time  resolved  anode  current.  Operating  con¬ 
dition:  Va  =  450  V,  rna  —  8.5  seem,  I*.  =  0.5  A, 
rhc  =  1.0  seem.  This  occured  while  the  thruster  was 
warming  up,  and  then  transitioned  to  the  high  cur¬ 
rent  mode,  as  in  Figure  16(e) 
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(c)  Time  resolved  anode  current.  Operating  condi¬ 
tion:  Va  =  450  V,  rha  =  6.5  seem,  I*.  ■  1.5  A, 

ill c  =  1.0  seem. 


0.00E+000 
Tine  (sec) 


(e)  Time  resolved  anode  current.  Operating  condi¬ 
tion:  Va  =  475.7  V,  rna  =  5.8  seem,  I*.  =  1.48  A, 
rri c  =  1.0  seem. 


Frequency  (Hz) 


(b)  Time  resolved  anode  current.  Operating  con¬ 
dition:  Va  =  450  V,  rha  =  8.5  seem,  I*.  =  0.5  A, 
rhc  =  1.0  seem. 


(d)  Time  resolved  anode  current.  Operating  con¬ 
dition:  Va  =  450  V,  rha  =  6.5  seem,  I*.  =  1.5  A, 
rhc  =  1.0  seem. 


Frequency  (Hz) 


(f)  Time  resolved  anode  current.  Operating  condi¬ 
tion:  Va  =  475.7  V,  rha  =  6.8  seem,  1^  =  1.48  A, 
rhc  =  1.0  seem. 


Figure  17.  Time  resolved  anode  current  for  the  low  current  mode. 
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IV.  Conclusions  &  Future  Work 


The  plume  measurements  expand  upon  and  confirm  previous  RPA  and  faraday  cup  data.6  A  comparison 
of  RPA  and  Faraday  data  for  identical  operating  conditions  shows  that  the  region  of  the  plume  with  the 
greatest  current  density  contains  the  most  energetically  uniform  ion  populations.  However,  the  RPA  results 
indicate  that  highly  energetic  ions  are  expelled  at  large  angles  relative  to  the  DCFT  axis  of  symmetry,  though 
the  current  densities  at  these  angles  are  much  lower.  Time  resolved  anode  current  data  reveal  a  distinct 
difference  between  the  low  and  high  current  modes.  The  data  demonstrate  that  the  anode  current  in  the  low 
current  mode  is  much  less  oscillatory  than  in  the  high  current  mode.  Identifying  the  cause  of  this  difference 
in  dynamic  behavior  will  likely  contribute  toward  the  determination  of  basic  DCFT  operating  principles  in 
each  mode. 

While  running  the  thruster  for  the  purpose  of  this  investigation,  the  operating  conditions  of  the  DCFT 
and  cathode  were  found  to  be  strongly  coupled,  motivating  a  future  study  that  systematically  examines  the 
effect  of  cathode  location  and  operating  condition  on  the  performance  of  the  DCFT.  A  detailed  look  at  the 
ion  energy  distributions  at  high  angles  in  the  low  current  mode  would  also  be  beneficial,  to  determine  if  high 
energy  ions  are  as  prevalent  at  high  angles  in  the  low  current  as  in  the  high  current  mode.  Additionally, 
repeating  some  of  the  measurements  made  hear  at  different  chamber  pressures  would  address  the  role  that 
charge  exchange  collisions  with  background  neutrals  play  in  deflecting  ions  to  high  angles.  However,  it  is 
likely  that  the  strong  magnetic  field  outside  the  thruster  is  a  strong  contributor  to  the  high  ion  energies 
observed  at  large  plume  angles. 
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